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Dr.  Edwin  P.  Rood 

Fluid  Dynamics  Program  (Qode  1132F) 

Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  VA  22217 

Subject:  Progress  R^xort  on  Contract  #N00014-91-J-4058 
"Ccnputational  Modeling  of  Marine  Prcpulsors'* 

For  Period  of  January  1,  1992  to  March  31,  1992 

Dear  Dr.  Rood: 

Ihis  letter  r^resents  our  third  qi:Brterly  progress  r^xort  on  the 
subject  contract  for  the  indicated  period.  In  addition  to  planned 
chjectives  stated  in  our  last  quarterly  r^xort,  namely,  cxjntinued 
efforts  on  the  3-D  parametric  Navier-Stokes  study  of  the  curved  blade 
passage  problem  and  further  validation  of  our  3-D  panel  cxxde  for 
marine  prcpellers  with  unsteady  non-uniform  inflow,  cur  focus  during 
this  period  included  extensions  of  our  Navier-Stoki^  code  to 
axisyinmetrical  flew  with  swirl  to  study  the  inflow  velocity  profile 
for  rotor-stator  blade  design  and  to  a  time-accajrate  Navier-Stokes 
solver  for  2-D  unsteady  flew.  We  have  also  initiated  a  numerical 
analysis  of  acoustic  noise  produced  by  2-D  blade-vortex  interaction. 

A  brief  summary  of  the  progress  in  these  areas  follows. 

2-D  TIME-ACCURATE  NAVIER-STOKES  SOLVER.  A  2-D  time-accurate  N-S 
solver  was  developed  using  artificial  cenpressibility,  (dual  time¬ 
stepping  and  dynamic  gridding  to  accemmodate  the  change  of  the 
ccsrputational  desnain  with  time.  Demonstrative  calculations  for  a 
pulsating  pipe  were  conducted  and  cenpared  with  available  perturbation 
solutions.  Ihe  configuration  of  the  pipe  and  the  pressure 
distributions  on  the  wall  are  presented  in  Figs,  la  and  lb. 

2-D  ACCUSne  NOTSF.  PRORTFM.  Based  on  OUT  2-D  unsteady  panel  cxxde,  an 
acxxistics  cuialysis  of  a  blade-vortex  interaction  flew  problem  (see 
Fig.  2)  was  initiated  using  an  cisynptotic  ea^ansion  techniejue  and 
small  perturbation  formulaticai.  In  this  example  problem,  a  vortex  is 
released  upstream  and  its  position  is  trached  as  it  is  cxxnvected  past 
an  airfoil.  During  this  time,  the  flowfield  is  calculated  at  each 
time  step  and  along  with  tracdcing  the  new  locations  of  the  initial 
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vortex,  the  wake  vortices  shed  fran  the  blade  trailing  edge  are  also; 

obtained.  Hiis  will  allow  xis  to  include  the  effects  of  the 

oontributiQn  to  the  acoustic  field  fraa  the  change  of  tlie  vortex 

distributicai  cti  the  blade  surface  as  well  as  that  due  to  the 

relocaticai  of  the  vpstream  vortices,  both  of  vhich  are  caused  by  the 

blade-vortex  interaction  vhen  the  blade  passes  by  the  vortices. 

Previous  calculaticxis  of  this  problem  have  ignored  the  induced  vortex 

effect.  Work  in  this  area  is  progressing  well  and  will  be  continued 

in  the  next  r^xirting  period. 


AXISyMMETtRICAL  FI£W  WITH  SWIRL.  Frcm  experimental  measurements  it  was 
found  that  the  inflcw  velocity  profile  for  fen  blade  design  considered 
earlier  is  not  uniform  in  the  radial  direction.  By  emitting  this 
radial  variaticxi,  the  total  flow  volume  can  be  overestimated  by  as 
much  as  15%.  It  was  further  confirmed  throu^  numericcil  anedysis  that 
this  radial  variaticxi  of  the  inflow  velocity  is  not  caused  by  viscous 
wall  effect  either  on  the  hub  or  on  the  casing.  Ihus,  it  was  hi^ily 
possible  that  the  non-unifozm  distribution  of  the  inflow  velocity  is 
caused  by  centrifugal  force  due  to  swirl.  A  N-S  solver  for 
axisynmetriced  flow  with  swirl  was  develcped  ard  code  validation  for 
siitple  geometry  was  conpleted.  Oenputations  for  realistic  geometries 
(see  Fig.  3)  are  currently  in  progress  and  will  continue  during  the 
coming  contract  period. 


3-D  PARAMETRIC  STUDY  OF  NAVIER-STOKES  FlfW.  As  used  in  previous 
studies  by  many  ether  people,  a  curved  rectangular  duct  was  chosen  as 
r^resentative  of  the  secondary  flows  in  turbcmachinery.  A  parametric 
study  of  the  curved  duct  problem,  including  aspect  ratio  of  the  cross 
section,  the  Deem  number  emd  the  Reynolds  number,  are  currently 
underway  and  their  results  are  yet  to  be  summarized.  At  present,  an 
algebraic  turbulence  model  is  used  in  our  3-D  N-S  code  for  hi(^ 
Reynolds  number  flow. 


PLANS  FOR  TOE  OCMING  PERIOD.  During  the  coming  three  months,  we 
anticipate  continuing  efforts  in  edl  the  above  areas.  Our  ertphasis 
will  be  on  ocnpleting  the  acoustic  einedysis  for  the  indicated  vortex 
problem,  and  starting  to  extend  to  more  ocnplicated  problems.  Work  in 
the  unsteady  N-S  aneilysis  for  2-D  airfoil  flow  problems  will  continue 
and  the  study  of  axisymmetric  flow  with  swirl  should  be  completed.  In 
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addition,  our  initial  catputations  of  the  3-D  parametric  study  of  N-S 
flow  will  be  summarized  and  a  low  level  effort  will  be  lo^  to  further 
validate  our  3-D  unsteady  panel  code  for  marine  propeller 
applications. 


Sincerely, 


CUVpm  Charles  L.  Merkle 

Distinguished  Alumni  Professor 
of  Mechanical  Engineering 
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Figure  3.  calculation  domain  for  fan  design 


